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ABSTRACT 


~~ 
Bee ea ‘ ‘ . P . . : 
~This investigation considers a non-linear technique for 


extracting natural resonances from transient electromagnetic 
scattering responses of radar targets. These natural reso- 
mances represent the complex poles of the target's transfer 
function in the Laplace transform s-plane. The advantage of 
their use in identification is their dependence only upon the 
geometry ard composition of the target and not upon the aspect 
and polarization of the incident signal. 

Based on recent theoretical research, a new signal model 
has been developed for describing the form of the transient 
scattering response. This form precludes the optimal use of 
previous methods for natural resonance extraction based on 
Prony's algorithm. In this effort, a modified least-squares 
approach is taken which can accommodate the actual form of 
the transient signal. The algorithm is tested using simulated 
signals with various noise levels and conclusions are drawn 


concerning the viability of the method. _ 
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I. INTRODUCTION 


This investigation considers a non-linear technique for 
extracting natural resonances from transient electromagnetic 
scattering responses of radar targets. These natural reso- 
Mances represent the complex poles of the target's transfer 
function in the Laplace transform s-plane. The advantage of 
their use in identification is their dependence only upon the 
geometry and composition of the target and not upon the aspect 
and polarization of the incident signal. 

Concepts and methodologies evolving from this idea have 
been developed based on the Singularity Expansion Method 
(SEM) by Baum [Ref. 1] in which the location of s-plane 
singularities characterizes a target. The direct extraction 


of poles (and residues) from the complex impulse response was 


initially attempted by Main and Moffat [Ref. 2]. These two 


“e 


early references spawned the idea of using natural resonances 
as the basis of non-cooperative target recognition (NCTR). 
The application of the Householder orthogonalization 


method using a successive orthogonalization process and the 
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eigenvalue method, that finds eigenvectors and eigenvalues, 


can determine the number and value of poles from the tradi- 
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tional transient signal model of a scatterer if there is a 


a 
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sufficiently low level of noise in the data [Ref. 3]. But 
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both methods show poor quality of performance if we apply 
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them to the new signal model, even were there is no noise 
assumed in the signal data. This new signal model, which 
represents the electromagnetic scattering transient response 
of actual targets, was recently derived by Morgan [Ref. 4]. 
As considered in [Ref. 2], the relationship between a 
scattering target and the singularities being extracted from 
the natural complex resonance signal has been emphasized in 
the context of the number of poles and the accuracy of poles 
to ensure the unique relationship between an object and its 
image. Namely, the accuracy of pole extraction plays an 
important part in setting baselines to check consistency 
between a received signal and a set of singularities. Also, 
(Ref. 2] describes the modeling problem of a signal such that 
we express a transient response signal in terms of complex 
exponentials and the way it can be manipulated by Prony's 
algorithm. This signal model (termed here as the traditional 
model), was found not to represent the actual transient signal 
in the recent work by Davenport (Ref. 5]. In 1982, Manilha, 
in his work, commented that the new signal model might be 
partitioned into “early time signal” and “late time signal", 
each of which characterize the signal by the time instant at 
which the traditional model can be acceptable or not [Ref. 3]. 
This effort considers the development of a method that is 
capable of handling the “early time" whose time varying 
residues of exponentials do not become constant until the 


excitation of the incident field disappears. Such a signal 


SE IRE SS aU a Rear eon ar et ee en 


can not be adequately or accurately modeled by the tradi- 
tional model. 

It is the intention of this work to investigate a method 
that can handle the “early time signal" and show "Robustness" 
under relatively heavy noise pollution as well as improving 
the accuracy of the poles we are int e.ted in evaluating. 

A non-linear parameter optimization wJproach, through the 
modified least-squares minimization - od, has been tried. 


In particular, we have been interested in evaluating its 


as oe on a eS ee ee 


performance in improving the accuracy of the poles extracted. 


acta 


The non-linear parameter optimization approach using itera- 


“sR Ra 


tions has contribued to improvement of accuracy and showed °- 
the possibility of application to the "early time signal" as 
was the final goal. We define the "new signal model" as a 
causal connection of the "early time signal" and "late time 
signal" throughout this paper. The new signal model will be 
presented in Chapter II, with both the description of the 


transient response mechanism and the way in which the signal 


model is constructed through the transient response mechanism. 


Chapter III briefly describes the problems associated with the 
traditional model, and with many of the current extraction 
methods in existence. Chapter IV is devoted to the descrip- 
tion of the non-linear parameter optimization through a 
modified least-squares method and the way in which the algo- 
rithm is implemented. Chapter V presents test results of an 


attempt to use this approach for pole extraction and its 


rare! he 


performance characteristics. Chapter VI contains a summary 
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a of the simulation results, the problems encountered, the 
potential problems expected, and the possibility of 
application of this method. 

Computations were performed using the IBM 3033 system at 
the Naval Postgraduate School. The testing of concept feasi- 
bility was the major impetus behind this effort, with concerns 
for practical implementation in a "real-time" environment left 
for subsequent study. Considerations to the processing speed 
will not be addressed in this thesis. All the calculations 


were done in a double precision environment. 
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II. SCATTERING TRANSIENT RESPONSE SYSTEM 


A. TRANSIENT RESPONSE MECHANISM 

An electromagnetic wave incient upon a scattering body 
forces currents to be distributed such that Maxwell's equa- 
tions and the corresponding boundary conditions are satisfied. 
These induced currents produce a scattered field that can be 
defined at any coordinate in free space if the expressions for 
those currents are made correctly. Changes in the angle of 
incidence, the incident wave shape, the polarization mode, 
and the target geometry and composition affect the results of 
the evaluation of the scattered transient response. 

Basically, the singularities uniquely characterize a 
scatterer by their number and locations. The corresponding 
residues, which are interpreted as the weights of these poles, 
either can be viewed in the frequency domain or in the time 
domain. 

As shown in [Ref. 3], equation (2.1) represents the form 
of the scattered signal impulse response of the target, in 
digital form 


H(k) = Hy (kT) + a: exp (kS_T) (2.1) 


where H(t) = § fort > Tg, T is the sampling interval and 
the infinite sum represents the natural resonance response 


with complex Se Oat jo, in the left-half plane. The poles 


m 
constitute the parameters for NCTR of an object. The process 
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first collects the data points within a finite time period 
(time window). 

The data collected is represented using equation (2.1) 
and equation (2.3) which partitions the signal into “early 


time” and “late time” components. 


oo M-1 
x(k) = | H(k) * J I(k) (2.2) 
k=0 k=0 
Tg-1 M=-1 © M~1 
X(k) = } H(k) * J I(k) + J] H(k) * J I(k) (2.3) 
k=0 k=0 k=T¢ k=g 
early time data late time data 


where I(k) is the incident wave at t = kT 


T is sampling interval in time. 


We can write the equation (2.3) in the form of equation (2.4) 
or (2.5). Both of these expressions display the essence of 
the new signal model which will be discussed in the next 


section. 


N 
X(k) = L An (k) exp (S_k) (2.4) 
n= 


where A, (k) = A, for k > Tf 


“ 


Naas N 
es X(k) = E(k) + } A exp(S_k) (2.5) 
ee n=l n n 

ey 

ey 

er = > 

‘ori where E(k) = @ for k > Tf 

ay and A, is constant for all k 

ae 
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T@ is the point at which the late time signal starts, so we 
consider that a time window might contain the data points 
ranging from the one before Tg to that after Tg. If we 
perform processing on a data window containing the data 
points which were from that early time region, we may lose 
some of the poles of natural resonance what are likely to be 
seen in the late time region. The data points after Kg 
presumably has all the information regarding poles, even 
though there are some extra poles that may be discriminated 
from the natural poles. These extra poles are introduced 
by the result of noise with relatively low SNR. Also, there 
is another risk to lose the natural poles having very high 
damping coefficients. There is no effective method yet found 
: to compromise on having the point in which we can see all the 
poles with less extra poles under the higher SNR condition. 
Meanwhile, the discrimination of E(k) appears to us as 
the quantity to be handled to give improvement in accuracy 
of the natural poles. The next section describes the new 
Signal model based on these aspects and the transient response 


mechanism. 


B. IMPULSE RESPONSE SIGNAL MODEL 


3 8 
Ph 


4, 
LACEY Yo, 


The late time impulse response can be expressed as a 


summation of exponentially damped sinusoids as in equation 


(2.6). 
ry bs N 
ae | X(k) = Y A, exp(S.k), k=9,T,2T,..., (M-1)T (2.6) 
oe, n=l 
ne: 
ae. 14 
A, 
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where T = At, and M is the number of sampling points. 


The signal model expressed in equation (2.6) is valid after 
the incident wave has completely illuminated a target and 
no forced response remains. 

As the mechanism of a transient response system is strictly 
governed by the physics and orientation of the target and the 
attributes of the incident wave, we can write the equation 
(2.7) as an implicit form of equation (2.2), assuming we 
strike the target with a plane wave impulse. 


© Tg-1 
x(k) = ) - 2 [ulk) - ulk-T9) J (2.7) 


Combining both equations (2.6) and (2.7), equation (2.9) would 


be sufficient if we assume N poles are present. 
X(k) = X, (k) + X, (k) 


0 dia 
X(k) = } H(k) ° [u(k) - u(k-T) ] 
k=9 k=9 


+ A, exp (Sik), k=$,T,2T, ..., (M-1)T (2.9) 


We call the impulse response signal model, in the form of 
equation (2.9), an implicit form of the new signal model. 
A simpler form of equation (2.9) can be either of equations 


(2.19) or (2.11). 


N 
X(k) = E(k) + L An exp (S_k) (2.19) 
n= 
N 
X(k) = } A, (k) exp (S_k) (2.11) 
n=l 
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for k > Tf 


where A, (k) = Ans 
In equation (2.1), the residues are time-varying until 
the time at which the augmented function E(k) in equation 
(2.19) becomes zero. Thereafter the A, (k) = A, (constant) . 
To begin with equation (2.19), it is necessary to define 
a new independent variable as a parameter whose time behavior 
is not predictable along the positive time axis. It is con- 
venient to define this as a set of variables. By regarding 


them as independent, the equation (2.11) would be in the form 


of the equation (2.12) in the discrete digital data processing 


sense e 
N 
X(k) = E(k) + } A, exp(S_k) (2.12) 
n=1 
T#-1 N 
X(k) = J e@.(k) + J A, exp (S_k) (2.13) 
n=] ® n=1 a 


assuming e,, (k) = Eo (k-n) 


The equation (2.14) is expressed in terms of a positive 
counting sequence ink. As in the equation, we might evaluate 
the equation (2.13) if we know the exact number of poles. But 
we know that it is impossible because we are in the region of 
the early time. So we write the equation (2.14) as a time 
shifted version such as the equation that has meaning from 


the observability viewpoint. 


ee 
X(k) = e_(k) + A. exp(S_k) (2.14) 
n=l ® n=] ® = 
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where k=T9+m, Td+m-1l, ..., TH+1,TH,TH-1, ... , TH,G 


for m > 0 
Here, the examples of the new signal model are presented 
in Appendix C, in the form of a decomposed signal. This syn- 
thetically generated signal is written in terms of e, (kK) 's 
and the sum of exponentials. In Chapter V, we present the 
decomposed signal from the synthetically generated data signal 
and that obtained by constructing the E(k) as the results 


of the non-linear parameter optimization processing. 


C. NOISE 

White Gaussian noise is assumed throughout this work. By 
introducing the noise into the new impulse signal model, we 
call the following equation the modeling function of a tran- 
sient response system and the complete and general form of 
the new signal model. 


SS N 
X(k) = N(k) + J e(k) + y A exp (S_k) (2.15) 
n=1 n=1 7 
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III. PROBLEMS WITH THE TRADITIONAL METHODS 


In this chapter, a brief description of problems with 
the classic Prony's method, the Householder orthogonalization 
technique and the Eigenvalue method are provided. Each of 
these methods computes the location of poles in the s-plane 
but only the last two methods can estimate the number of poles 


prior to the discrimination of the poles. 


A. PRONY's ALGORITHM 
Equation (2.3) is interpreted as the following difference 
equation assigning the values of the real part of the signal at 


specified sampling points in its left side of the expression. 


(3.1) 


N 
X((M-1)k) = } AQZ) (M-1)k 
n=1 


And a polynomial equation (3.2) that has the same roots 2, 


N n 
ata a,Z2 = 6 (3.2) 


18 


can be combined with equation (3.1) to yield the Prony's 


difference equation in the form of equation (3.3) 


N 
Y. a,X(n+k) = 0, for k=0,T,2T, ... , (M-1) T (3.3) 
n=l 


Then, we write equation (3.4) from equation (3.3) as 


N N-1 
dank (atk) = Lo a,X(ntk) + ayx(n+k) = 9 (3.4) 
n= n= 


Using matrix notation, the equation (3.3) may be written as 
equation (3.5) or (3.6). 


N-1 
a a,X(ntk) = -a,X(n+k) (3.5) 
nz 


Let ay = 1, then equation (3.5) becomes as 


Xy-1 Ay.1 = *y (3.6) 
where Xye1 is N by N circulant matrix of sampled data 
Ay.1 is N by 1 Prony's coefficients matrix 


Xy is N by l row matrix of sampled data set 


[Xgenene coe * Xan] 


The Prony's coefficient a,'s are to be calculated by 


Making use of the characteristics of the circulant matrix X. 


ay 
Ay-1 = (Ayia Xye.) ~ Ayaan ¥n (3.7) 


As described in the above, this method can be applied after 
the number of poles are known to us, then the value of the 


s,‘s are to be found from the equation (3.8). 
19 


S,, = LN(Z,)/k, for k > 2N (3.8) 


B. HOUSEHOLDER ORTHOGONALIZATION METHOD 


\ 
> A general expression of equation (3.2) may be represented 
os by equation (3.9) assuming the number of poles are not known. 
1 Ni 
ag + a2 Hoss SO an = g (3.9) 
where N‘' is the unknown variable 
Then, the equation (3.9) is to be satisfied under the condi- 
tion that there are enough sampled data points m > 2N. 
agX, + a,x) +... t ayy = g (3.19) 
where the N by 1 sampled data matrix is defined as 
Application of the successive orthogonalization through the 
Gram-Schmit process would produce the orthogonal vector set. 
0 = fog, Ors see On} 
n-1 
where 0, = x, ~ BE) <K10,> 0; (3.11) 
i=9g 
If we set the Og to be 1, then the N by 1 Xn Matrix is going 
to be in the form of equation (3.11). 
-1 
a) z ” 
- Hy 7 On + 1 Meye Oy > Oy : soon) 
hove, i=g 
es 
Peet So that the equation (3.12) holds as X=0M. In the above, 
Ns 
we M is the 2-D matrix of multiplication factors whose diagonal 
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elements are all l's. Here, an orthogonal vector oO, makes 
the corresponding data set x, to be orthogonal to all the 


previous vectors xg - x If we have had all the x vectors 


i-1° 
during i-l times of successive orthogonalization process 
such that there is a non-zero orthogonal vector. We have to 


receive continuously the next data set x, and check whether 


i+l 
the orthogonal vector in the next step is zero or under a given 
threshold. If it vanishes, we can say that we have (i-1) 

poles. Our test against noise polluted data signal using the 


Householder method showed the inability of handling the early 


ewe we eee - 


time signal as well as the data being heavily polluted. 


C. EIGENVALUE METHOD 


Again, from the equation (3.19), we can rewrite this 


. equation as in the form of equation (3.14). 
Xyr2ye = g (3.14) 


where x = [xg] x, | Carte 1X0] 


The eigenvector can be derived from the equation (3.15) by 


a* 
Ba 


=. 
a 


mt 


doing some mathematical manipulations of equation (3.14). 


Gn: ary ae = i ee = § 


bat Wass Sake are 
544% 4,4, Ee 


d From this equation, we find the eigenvalues correspond to 
a eigenvectors to see if there is any eigenvalue approaching 
zero. If there is one zero or under threshold, we also 


f. 
q consider the number of poles as N'-l. Even both the 
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Householder orthogonalization method and the Eigenvalue 


method can provide the means of calculating the number of 
poles, their algorithm (tending to fit the data points to 


the traditional signal model) showed the lack of generality. 


IV. NON-LINEAR PARAMETER OPTIMIZATION APPROACH 


A. INTRODUCTION 

In 1963, Marquardt suggested an algorithm for the least- 
squares estimation of non-linear parameters [Ref. 6] when 
highly accurate parameter values are required. As we can 
see in the signal model in equation (2.14), that new signal 
model has multiple parameters that are functions of time. 
Here, we rewrite that equation again in the form of equation 


(4.1). 


Tg-1 N 
X(k) = N(k) + ; e (k) + 5 Ay exP(S)k) (4.1) 


n=k eS 


n=1 


° 


where e,, (k) = 9, for k >T% 
We are not interested in the value of .@, (k) itself, but 
include their estimations in order to contribute to the 


accuracy of s,'s and A'S: 


In other words, the extrapolation of the sum of exponen-~ 


| 
bat 


tials, which is assumed to have all the poles, to points in 


the early time region may provide more accurate poles and 


TRB: 


residues simultaneously. 


2 
WA 


The major advantage is the fact that we can make use of 


ars? 


hans 


the high SNR data signal in the early time region with the 


basic information of poles that were derived in the late time 


ats ud 


é: 
a 


region. 
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We will be processing on a data window whose vector 
length is to be increased one by one by moving the first 
element of that vector toward the time-origin point. The 
optimization process evaluates the normal equation to find 
the local optimized set and finds the global optimized set 
as its final goal. 

Let us define an ERROR function as in equation (4.3) 

N-1 

ERROR{A,,S,,E(k)} = y [X'(k) - x(k) ]? (4.3) 

n=1 

where X'(k) is a measured data point 
To have minimized the least-squares error at every instant 
of measured time, the ERROR function has to be minimized by 
obtaining the global set of parameters being optimized. We 


also can write equation (4.3) in a more concrete manner as in 


equation (4.4). 


M-1 2 M=-1 2 

ERROR {xX} = J [X'(k) - X(k)]°= } e(k)e (4.4) 
k=9 k=9 

where X = [E(k) ,A,,S.] 


Non~linear parameters such as e, bk) must be obtimized in 
a way such that the more accurate values of S. have to be 
calculated as we increase the number of non-linear parameters 
e,- Users may define the numbers of parameter e,'s that 
are at least equal to or greater than the number of data 
points in the early time region. Now it is our task to find 
out the optimized value of AneSne and en through optimiza- 


tion processing, either in the global sense as it is represented 


in the equation (4.5) or using the normal equation (4.6). 
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9 ERROR(X) /dX 


it 
we 


(4.5) 


0 
Qa 


de, (k) 8 X(k) (4.6) 


where k=9,T,2T ... , (M-1)K 


n=1,2, ..., N 


With using the normal equation, as it is shown in the (4.6), 
a global set of non-linear set of parameters can be obtained 
through an iterative evaluation and fitting process. In our 
simulation work, the modified least-squares method was the 


basis of the equation (4.6). 


B. THE MODIFIED LEAST-SQUARES METHOD 
The error function is rewritten as in the expression of 
the equation (4.7). 
ERROR{K,x} = X'(k) - X(K, xX) (4.7) 
where X'(k) is the k-th sampled data 
X(k, x) is the model function 


xK is a vector containing the k-th sampled 


V] 


Let us define xn to be an initial estimated value of x, then 


a sequence of approximations to the optimized value is to be 


generated by the equation (4.8). 


cae m+l m 


cs T -1l .T m 
e: x =x - [a,D., + Jn Jnl Jin ERROR(x ) (4.8) 
. one 
tw where Tn is the numerical Jacobian matrix evaluated 
ats 
“te m is iteration number of successive optimization 
or 


Da is a diagonal matrix equal to the diagonal of 
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an is a Marquardt parameter 


The number of total iterations can be controlled by the 


threshold which is defined as in the equation (4.9). 


m+1 m (4.9) 
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V. TEST RESULTS AND PERFORMANCE EVALUATION 


A. INTRODUCTION 

In order to establish the ability of the program listed 
in Appendix B to extract the poles and residues correctly, 
three different simulated signals, each of which is polluted 
with infinite, 30 dB and 15 dB SNR, were created by the syn- 
thetic signal data generation routine. These 3 sets of signal 
data were chosen to span a range of possible situations. 

In the context of the transient signals and additive sta- 
tionary Gaussian noise being used, the definition of SNR is 
in terms of a ratio of energy quantities intergrated over the 


entire 20 nsec time window. 


TABLE I 


SIMULATED SIGNAL 1 


POLES 
RESIDUES (Nep. GH2) 
1.9 + il .@ -1.9 + i1.9 
1.9 - il.g -1.9 - i1.9 


Simulated signal 2 consists of 2 sets of pairs of complex 
conjugate poles and residues, which is an extrapolation of the 
simulated signal 1, using its parameters. Simulated signal 3 
has 3 sets of pole-residue pairs, extrapolating from signal 2. 

There are 5 options in the signal generation program in 
choosing a particular function of E(k). In this simulation 


work, a trapezoidal wave form of E(k) was used with T#=9 .42 


nsec. So that the unknown parameters that are to be augmented 
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at every processing step will be a maximum of 10 plus 4 times 
the number of poles at the final processing stage, when we 


have 512 data points within a 20 nsec time window. 


TABLE II 


SIMULATED SIGNAL 2 


POLES 

RESIDUES (Nep. GHZ) 

1.9 + i1.6 -1.9 + il.g 

1.9 - i1.g -1.9 - i1.g 

g.5 + ig.5 -2.8 + i2.g 

9.5 + if.5 -2.8 - j2.9 

TABLE III . 
SIMULATED SIGNAL 3 
POLES 
RESIDUES (Nep. GHZ) 

1.9 + i1.f -1.9 + i1.g 
he 1.g - il.g -1.g - il.g 
sy §.5 + ig.s -2.9 + i2.9 
“ao 9.5 - ig.5 -2.f - i2.9 
ee 0.25 + i9.25 -3.9 + i3.9 
- $9.25 - ig.25 -3.8 - i3.9 
ae The results of varying the number of additional data 
vA points from zero (all the data points are from late time 
Nt region) to ten extra data points (10 extra points are from 
me early time region are added to the points of late time 
a region) are contained in Table IV through XII. These nine 
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cases correspond to the 3 synthetic signals, each having 


3 additive noise levels. The 3 synthetic signals are plotted 


in Appendix C. In the next section, the accuracy of pole ’ 
extractions, as indicated in the tables are displayed by way 
of graphical pole maps. In addition, reconstructed waveform 
obtained from the parameter extractions are compared graphi- 


cally to the original waveforms. 
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TABLE IV 


PARAMETER OPTIMIZATION FOR SIGNAL 1 (NOISE FREE) 


Bee ee I eg Se a mee 


RESULTS 

TARGET TYPE sTGT-1 

wAVEFCRM TYFE SE CNFTR 

CONTACT CATE:CEC 15 

FILE NUMESFILEQO2 

NUMBe CF FOLE: 2 

TAGLE CF RESICUES AND POLES 
S3EB2SS 2252 SSS SSBZESES SEES SZSESERE 

PAIR # RES@-KEA RES~IMAG POLE=REAL FOLE=IMAG 
1 1.Q9¢¢€00000 1600000000 -1.000C0C000 1200000900 
2 1.CCCCIGCD «=-le CCOCCODO )=~-1.000CC00G -1.00000000 

RESGLTS CF OFTIMIZATICN wiTh O EXTRA CATA PCINTS 

PAIR # RES e*REAL RESe=IMAG POLE=-REAL FOLE-IMAG 
1 Ceo $S$S96610 1.00001592 -0.99994705 C-99$989 71 
2 009595596563 =-1600001575 -00999S4777 -0.99598971 

RESULTS CF GFTIMIZATICN WITH 2 EXTRA DATA PCINTS 

PAIK # RESeREAL RESe- IMAG FCLE=REAL FOLE~IMAG 
i 1-0CC00251 00$999S667 -1.G00CCO25 10C0000 85 
2 1.CCC902Z51 9 -0.59999667 -1.00000026 -1.00000085 

RESULTS CF OFTIMIZATICN WITH 4 EXTRA DATA PCINTS 

PAIR # RES «REAL RESe 1 MAG POLE~REAL FOLE=IMAG 
i 1.9CC013C8 1.0C900329 -1.000C153€ 0 $959995S 
2 LeOCCCI312 «==1eC00000330 -1200001545 -1400000092 

RESULTS CF OFTIMIZATICN wWITh 6 EXTRA OATA PCINTS 

PAIR @# RES es R55 RESe~IMAG POLE-REAL FOLE-IMAG 
1 0.566$285213 1600303040 -0259975425 0099S93208 
2 0696685210 ~1eCCO03047 -0.59975441 —-0.99998213 

RESULTS CF OFTIMIZATICN &ITH 8 EXTRA DATA PCINTS 

PAIR # RES .-REAL RESe~ IMAG PCGLE-REAL FOLE-IMAG 
i i. CC€C0036 1200001011 -1-00000472 099999818 
2 LeUCCIO026 -1e 00001041 -1-COOCCS540 -0259599838 

RESULTS CF OPTIMIZATICN WITH 10 EXTRA DATA PCINTS 

PAIR # RES esREAL RESe-IMAG POLE-REAL FCLE-IMAG 
1 1-Q0CCCO4E7 Ce $$9S99062 -06S995S674 1.2C00001 65 

1.00030453 3 -0.99999067 -0.9999S735 -1.00000168 
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TABLE V 


PARAMETER OPTIMIZATION FOR SIGNAL 1 
(SNR=30 dB) 


TARGET TYPES TGT=1 
wAVEFGRM TYPE :ECA20TR 
CONTACT CATE:CEC 15 
FILE NUME:FILEQO3 
NUMB. CF POLE: 2 
TABLE CF RESIOUES ANO POLES 
BBZBSESBSZISS SIESSSESSSZEVBWESZT AIST I= 
PAIR # RES=REAL RES=IMAG POLE-REAL FOLE=ILMAG 
1 1e9CCCO0CD 1600000090 -1.050CG0Id 1499009000 
2 1.d€CCOOCS =-1.CCOCCCOO =-i.cOCcCOOC § -1.00000000 


RESULTS CF CFIIMIZATICN wITH O EXTRA DATA PCINTS 


PAIR # RES esREAL FES e-- IMAG PQLE=REAL FOLE-IMAG 
1 1e01CCES572 0092546346 -0077E7649S2 14200370076 
2 1203315188 ~=1.01073823€ -1630247522 C.99197436 


RESULTS CF OPTIMIZATICN wITH 2 EXTRA DATA PCINTS 


PAIR # RESe=FEAL RESe~I MAG POLE=REAL FOLE-IMAG 
1 1.0107C336é 1eCC3E0284 -14C095116i Ce$9S793 04 
2 1.0107€235 -1.00360275 <-1.00951163 059979306 


RESULTS CF OPTIMIZATICN wITH 4% EXTRA DATA PCINTS 

PAIR #  RESe-KEA RESe =I MAG POLE=REAL FOLE-IMAG 
1 120106€21 1.00307€55 -1.CQ095284S 0.$S9S77751 
2 1-010€€E32 ~1200367409 -1.00953782 3~-0299977735 


RESULTS CF OPTIPIZATICN wWITh 6 EXTRA DATA PCINTS 


PAIR #  RESeFEAL RE Se IMAG POLE~REAL F 
1 1.0108S740 1200357123 -1.0093C0702 Oo 
2 1.010£¢679 - 1200257126 -1.C0930751 Ce 


RESULTS CF GFTIMIZATICN WITh & EXTRA DATA PCINTS 


PAIR # RES es Ecal FES o- I MAG POLE@RcA FOLE-IMAG 
1 1.6106764 1.00 362571 ~1.0095058 0.999786 86 
2 120106€C35 - 1.00 262636 “Le COSECSEE =-06.99S7869S 


“. RESULTS CF QPTIMIZATICN WITH 10 EXTRA DATA PCINTS 


Pata # RESSTEEAY GES scImag  -FCLETAEAL FOLECIMAG 
a 1 1.010€¢€Z1 1.002€0375 =“1eCOSECSTE 00599579175 
2 2  201cetedd -t og seozes =1.6095063C —-04S9979177 
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TABLE VI 


PARAMETER OPTIMIZATION FOR SIGNAL 1 
(SNR=15 dB) 


TARGET  TYPE:TGT=1 
WAVEFGRM TYPE :HONISTR 
CONTACT CATE:CEC 20 
FILE NUME SFILEOO4 
NUMB. OF POLE: 2 
TABLE CF RESIOUES AND POLES 
BSESSBSSIEZ SZSASBASISESHRBSSTTZEBZASTA.Z 
PAIR # RES-REAL RES— IMAG PGLE-REAL —_ FOLE-IMAG 
1 leCCCSGOCO §«=— 1 €6000090 ~-1.00000005 1400009000 
2 1.0€CO0000 -1.C0000000 -1.00cCCO0C -120Ga00000 
RESULTS CF OPTIMIZATICN WITH O EXTRA DATA PCINTS 


Ses CEP BP 222 SSB QoeeQe 2eEBQ S22 22222 BeBe BWeaenwe 


PAIR # RESe REAL RESe~I MAG PCLE=REAL F 
0.5€C47SE2 Le 26922892 ~-02578SC08% Oo 
2 1259727710 -0.80336947 -1.57012325 <-0. 


RESULTS CF GFTIMIZATICN WITH 2 EXTRA OATA PCINTS 


PAIR # RES REAL RES. ~I MAG FOLE-REAL FOLE~IMAG 
1 1666657059 14602070489 -1.055CES4S 0459869875 
2 1.0€€57$23 -16C62071795 -120530627€ -0.59869876 
RESULTS GF OFTIMIZATICN WITH 4 EXTRA DATA PCINTS 
PAIR # RES REAL RESe-1MAG POLE“-REAL _POLE=IMAG 
i 1.06€5338 1.02072970 —~-1.05304404 0299868990 
2 1eO€CE2818 -1.02073354 -14052C40S1 ~0.99869088 
RESULTS CF GETIMIZATICN WITH 6 EXTRA DATA PCINTS 
PAIR # RES e=REAL RES. IMAG POLE“REAL _POLE-IMAG 
“y LeOECS1459 3602673547 -1.0525526€ 0699673653 
2 130€054674 <-120207279% 10053213212 ~CeSS864415 
RESULTS CF OFTIMIZATICN WITH 8 EXTRA DATA PCINTS 
nes PAIR # RES e-REAL RESe-IMAG PCLE-REAL | FOLE-IMAG 
“AN 1 1.02449 720 1. C¢88S916 0066485231 1201775335 
x 2 1207C48258 -1.10030318 -12039C7420 ~0.S30041 24 
a RESULTS CF OPTIMIZATICN WITH 10 EXTRA DATA PGINTS 
a went Ow SF 2 SESW G22 Bee 22 22s ea ween SB eaeewe 
PAIR # RESe-REAL RESe-1MAG POLE-REAL _FOLE-IMAG 
eae 1 10021665926 1611984550 -1.56675312 6.92761091 
a 2 10022£1816 -1507048377  -0065299821 -1202160947 
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TABLE VII 
TABLE CF RESICUES AND POLES 


SBSSBS 224 335 SESH SSERSREZAZSAE 


PARAMETER OPTIMIZATION FOR SIGNAL 2 
SULTS CF GCPIGMIZATICN hITH O EXTRA DATA PCINTS 
RES eo REAL 


- 
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RESULTS CF OFTIMIZATICN wITh 2 EXTRA CATA PCINTS 
RESULTS CF OFTIMIZATICN hITH 4% EXTRA DATA PCINTS 
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RESULTS CF OPTIMIZATICN &ITH 6 EXTRA CATA PCINTS 
PAIR # 


“lel 


RESULTS CF GFIIVIZATICN wITH 8 EXTRA DATA PCINTS 
RESULTS CF GPTIMIZATICN WITH 10 EXTRA DATA PCINTS 
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(SNR=30 dB) 


TABLE YIITI 
TABLE CF RESIQUES AND POLES 


PARAMETER OPTIMIZATION FOR SIGNAL 2 


RESULTS CF OFTIMIZATICN WITF O EXTRA DATA PCINTS 
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RESULTS CF OFTIMIZATICN WITH 2 EXTRA DATA PCINTS 


0e $9545224 


RESULTS CF OFTIMIZATIGN WITh 4 EXTRA OATA PCINTS 
RESULTS CF OFTIMIZATICN WITh 6 EXTRA DATA PCINTS 
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PAIR # 


8 EXTRA DATA PCINTS 
FOLE-REAL 
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RESULTS CF OFTIMIZATICN wITh 10 EXTRA DATA PCINTS 


RESULTS CF GFTIMIZATICN wITh 


RESICUES AND POLES 


TABLE IX 
(SNR=15 d3) 
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PARAMETER OPTIMIZATION FOR SIGNAL 2 


ae dO ee 


Ma Tate 


Tae 


8 Fete" 2 Tere 


DATA PCIN7T 


ANNES 


RESULTS CF OFTIMIZATICN wITh O EXTRA 
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RESULTS CF GFTIMIZATICN wWiTh 2 EXTRA DATA PCINTS 


POLE-REAL 


RESe-IMAG 
35 


RESULTS CF OFTIMIZATICN wITh 4 EXTRA DATA PCINTS 
Oe $3 

RESULTS CF OPTIMIZATICN wITH 10 EXTRA DATA PCINTS 
Ce $7? 


RESULTS CF OPTIMIZATICN wITh 6 EXTRA DATA PCINTS 
RESULTS CF OPTIMIZATICN wiITh 8 EXTRA DATA PCINTS 
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RESULTS CF CFTIVIZATICN wkITH 6 EXTRA DATA PCINTS 
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Figure 5.1. Estimated E(k) for Signal 1 of 15 dB SNR. 
(32 data points from late time region). 
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Figure 5.2. Estimated E(k) for Signal l of 15 dB SNR. 
(48 data points from late time region). 
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C. PERFORMANCE EVALUATION 

Unlike the traditional methods, this method revealed 
itself as one that can be applicable to a general scattering 
problem, including the early time data. A potential superiori- 
ty of this method over the traditional methods is found in its 
capability of handling the early time signal that has relative- 
ly large energy content and less sensitivity to the fixed noise 
level. 

Regarding the capability of handling the early time signal, 
the performance of this method seems to be very high compared 
with the traditional methods which do-not model that signal 
generically. But, two points need to be considered when we 
ensure this high degree of performance. One important point 
is that we be able to estimate properly the number of poles to 
optimally ask for. This will depend upon the SNR of the data, 
increasing generally with increasing SNR. The actual physical 
scattering data has an infinite number of poles but, because 
the energy content of natural modes tends to drop off rapidly 
beyond the bandwith of the excitation signal the higher order 
poles are ignored, blanding with the noise present. Another 
point is that the initial data window has to contain only the 
late time signal data. 

As the results of the tables in the previous pages indi- 
cate, converging characteristics of the value of parameters 
is not blocked by a high level of noise pollution. It has 


been observed that the computed value of poles are likely to 
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converge closely to the true values that are listed in Tables 


I through III either by choosing the number of poles as small 
or by increasing the size of the data window. Even though it 
is hard to prove this optimality of this method (due to the 
non-linearity of the scattering transient response signal 
model), handling of the early time signal seems to be cast 
into a simpler problem space and that optimality may be ex- 
plained by observing the parameters at every processing step. 
Namely, the parameter estimation of this method is less 
influenced by the unreliable points than the more accurate 
ones. Figure 5.1 and 5.2 and the results in the previous 
section explain the reason why the above assertion is valid. 
We see that the estimated value of parameters through optimi- 
zation is not going to worse as we include the early time 
data points (these data points can be regarded as the unrelia-~ 
ble data when we extrapolate the sum of exponentials of the 
late time region into those of the early time region) into the 
data window what we see. But, as we are adding more data 
points of the late time region, improved values of parameters 
are obtained through processing, as in the results and Fig. 
5.2. Although this is true, we still see that the plot of 
E(k) in Fig. 5.2 does not provide sufficient data to confirm 
that an improvement of accuracy depends only on the number of 
late time data points, because those data points have a lower 
level of SNR locally. Now, let us get back to Tables IV, VII, 


and X. For the signal inputs whose SNRs are assumed to be 
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infinite (noise free), we have a set of parameters whose 


values are not converging to the true values as we add an 
increased number of unknowns. These results ensure us that 
the increased computer round-off errors obtained by more 
data processing (more parameters to extract) will produce 
divergence. Even initially with low SNR, the increased 
signal energy in the additive early-time data may improve 
the results until, by adding enough extra data points and 
parameters, round-off.error will eventually catch up with 
the most optimized one and beyond a specific point the 
results will begin to diverge. So that, in a general sense, 
we can say that the improvement in accuracy may be obtained 
either by the number of data points or by making use of the 
data points of high SNRs. If we consider a signal having ; 
high damping coefficients we can not use the signal data 
what are likely to be hidden by noise. The technique used 
here in processing the data proceeds by sliding the first 
data point along the time-axis towards the origin until the 


data window expands its range up to time-origin by adding 


the unknown parameters E(k)'s which is the same as the number 


of data points in the early time region. 


Figures 5.3 through 5.11 shows how rapidly the parameters 
we want to derive converge to the true values, within a con- 


vergence bounds. Here, all the plots were reconstructed from 


the results. 


VI. SUMMARY 


A modified least-squares minimization method has been 
developed and tested against synthetically generated signal 
data using a non-linear parameter optimization algorithm. 

The effect of varying the number of data points either in 

the region of early time or in the late time under various 
typical noise environments has been studied, defining the 
Criteria for using this algorithm for the extraction of poles 
and residues under relatively heavy background noise, as well 
as for the time varying residue signals modeled by Morgan 
{[Ref. 4]. 

In an attempt to use this method for signals with high 
damping coefficients, 3 synthetically generated signals were 
used as the inputs (Appendix C). Extraction of poles from 
the direct synthetic signal waveform was done through non- 
linear numerical evaluations, not by the Prony-based evalua- 
tions which were used in the traditional method. The 
parameter estimating algorithm for this non-linear signal 
model had worked successfully, providing users with optimized 
parameters that were converging to the true values, within 
specified convergence bounds. It was found that although we 
have the early time signal data having very high SNR, those 
data can not contribute to improving accuracy by simply 
increasing the number of extra data points and extra unknown 


parameters. 


In order to get more accurate values of poles, we have 
to have an increased base of late time signal data, which 
have relatively low SNR. 

It has been shown that this method is "robust" even under 
heavy noise conditions. But three basic requirements have to 
be met for using this method as a general methodology; the 
optimal number of requested poles has to be known in advance 
to processing, reasonable initial estimation of parameters 
has to be made, and the transition time from early to late 


time signal models has to be known a priori. 
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APPENDIX A 
SYNTHETIC DATA GENERATION PROGRAM 
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NON-LINEAR PARAMETER OPTIMIZATION PROGRAM 
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